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General considerations in general relativity and quantum mechanics are known to potentially rule out contin¬ 
uous global symmetries in the context of any consistent theory of quantum gravity. Assuming the validity of 
such considerations, we derive stringent bounds from gamma-ray, X-ray, cosmic-ray, neutrino, and CMB data 
on models that invoke global symmetries to stabilize the dark matter particle. We compute up-to-date, robust 
model-independent limits on the dark matter lifetime for a variety of Planck-scale suppressed dimension-five 
effective operators. We then specialize our analysis and apply our bounds to specific models including the Two- 
Higgs-Doublet, Left-Right, Singlet Fermionic, Zee-Babu, 3-3-1 and Radiative See-Saw models. Assuming that 
(i) global symmetries are broken at the Planck scale, that (ii) the non-renormalizable operators mediating dark 
matter decay have 0(1) couplings, that (iii) the dark matter is a singlet field, and that (iv) the dark matter density 
distribution is well described by a NFW profile, we are able to rule out fermionic, vector, and scalar dark matter 
candidates across a broad mass range (keV-TeV), including the WIMP regime. 


I. INTRODUCTION 

Particle physics models achieve stability for dark matter 
(DM) particle candidates by advocating the presence of either 
discrete or continuous global symmetries. Discrete symme¬ 
tries arise, for example, from broken gauge (local) symme¬ 
tries, which are respected at the Planck scale lEia. Continu¬ 
ous global symmetries, instead, are generically violated at the 
Planck scale, leading to important implications on the dark 
matter phenomenology of the associated models. 

There are several reasons why continuous global symme¬ 
tries are not expected to be present in a consistent theory of 
quantum gravity, which rely on general facts in gravity and 
quantum mechanics: 

(i) No-Hair Theorem: Since local U{1) symmetries are 
effectively identical to Gauss’s law, any observer outside a 
Black Hole (BH) horizon can determine the BH charge. How¬ 
ever, if there existed continuous global symmetries, when a 
charged particle gets trapped inside the BH there would be no 
way to assess this from outside the horizon. Thus the charge 
would appear to be “deleted”, in contradiction to its conserva¬ 
tion ||3]. 

(ii) Hawking Radiation: The main problem with continu¬ 

ous global symmetries has to do with Hawking radiation H. 
Since there are no gauge interactions associated with contin¬ 
uous global symmetries, one could throw a large amount of 
charged particles into a BH and increase its charge (Q) in¬ 
definitely ilHlsl. However, the theory of Hawking radiation 
indicates that until Tnawking > tn, where m is the mass of the 
lightest charged particle pair, the BH does not radiate charge. 
Combining this with the bound on the BH mass, namely 
Qm < Mbh < Mpi/m, where Mpi = 1.22 x lO^^GeU is 
the Planck mass, we find Q < This limit can how¬ 

ever be violated by making Q sufficiently large. Hence, if Q 
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were conserved we could have identical BHs with an infinite 
number of states labelled by Q ^ 

(iii) Entropy: Since an external observer cannot infer a 
global charge, in order to assign an entropy to a given BH one 
would have to count all micro-states of all charges, finding 
an entropy of order ~ log(Q). Now, taking Q indefinitely 
large, one would violate the Bekenstein-Hawking formula, 
which says that entropy counts the number of states of a 
BH. Therefore, such objects are ruled out, as are continuous 
global symmetries 0. 

While there are general arguments for the breaking of 
continuous global symmetries at the scales of quantum 
gravity, those are not well-established. For example, Ref Q 
disputes such asguments based on non-thermal deviations 
from Hawking radiation. This notwithstanding , hereafter, 
we assume that continuous global symmetries are indeed 
broken at the Planck scale and show that the notion of such 
symmetries being broken at the Planck scale has profound 
implications on DM phenomenology. 

In this study, we assess the possibility of using continuous 
global symmetries to stabilize DM particles. In order to derive 
results applicable to a variety of particle physics models, we 
consider Planck-scale suppressed, dimension-five effective 
operators that mediate the decay of generic DM particles of 
spin 0, 1/2 and 1; the operators under consideration violate 
continuous global symmetries, and thus induce the decay of 
DM particles whose stability relies on such continuous global 
symmetries. The decay of long-lived but metastable DM par¬ 
ticles can inform us on the DM particle nature (see e.g. ClIHl); 
for example, stringent bounds on the lifetime of electroweak- 
scale DM stem from the observed diffuse gamma-ray flux 
0, which implies lifetime r > 10^® s, thus a billion times 
longer than the age of the Universe. We emphasize the fact 
that even though continuous global symmetries might break 
down to Zm discrete symmetries at low energies, one can 
always generically construct Planck-suppressed effective 
operators that would induce the decay of the DM particle: our 
results can thus be applied to any continuous global symmetry. 
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Naively, one might expect that Planck-scale suppression 
might have a negligible impact on the phenomenology of 
models which advocate the existence of continuous global 
symmetries to stabilize the DM particle. Using current 
cosmic-ray. X-ray, gamma-ray, neutrino and CMB data, span¬ 
ning the entire keV-TeV energy range, we show that, some¬ 
what surprisingly, continuous global symmetries are strongly 
disfavored as a mechanism to stabilize DM particles. In par¬ 
ticular, we rule out, under the aforementioned assumptions, a 
rather large DM mass range, including the classic WIMP mass 
range around the electroweak scale. 


II. OBSERVATIONAL CONSTRAINTS 

In this section we summarize how we derive our model- 
independent limits on the DM lifetime. We employ through¬ 
out our analysis an NFW profile iflOl . 

Pir) = , 

r/rs(l -I- r/rsy 

with a scale radius, Vg = 24.42 kpc, and ps = 0.184, such that 
to reproduce a local density of pQ = 0.4 GeV/cm^ ifTTI . We 
point out that our results would all scale linearly with other 
choices for the local dark matter density, unlike for the dark 
matter pair-annihilation case. 

I suspect the limits will essentially stay the same, since DM 
decay signals at high latitudes are not very sensitive to this. 


where H{z) is the Hubble rate and nniz) = 
is the number density of hydrogen in the Universe at a given 
redshift with = 1-9 x 10“^cm“^ being the present-day 
value, and 
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where /^ec is in general a function of the DM mass and 
redshift. We closely follow Ref. 1191 and average over the 
redshift dependence to get fdec as a function of the DM mass 
only, according to Table II of Ref ll20ll . We then modify the 
energy deposited into the intergalactic medium by inputting 
the equation above into the CosmoRec m package to 
compute deviations on the ionization history, which now 
depend on the DM mass and lifetime for a given decay final 
state. We obtained values for fd^c for the e'''e“,/i+/r”,T+T~ 
final states in agreement with Ref. Il20l . In addition, we 
calculated the efficiencies for decays into quarks and gauge 
bosons using the recipe of m. In detail, in order to compute 
the bounds for WW and qq final states, we rescale the limits 
from Ref |[T9l by using the appropriate energy deposition 
efficiency, which in those cases is slightly larger than for the 
pp final state, resulting into a slightly more restrictive limit 
as what shown in Fig. 1, and in agreement with ca. 


A. CMB data 

Precise measurements of the Cosmic microwave back¬ 
ground (CMB) provide robust limits on DM decays, since 
the latter alter the ionization and heating history of the CMB 
as well as its power spectrum. Using combined data from 
Planck im, WMAP9 ifTJI . Atacama Telescope llT4l . South 
Pole Telescope ca, Hubble Space Telescope IfThl and Bary- 
onic Acoustic Oscillations ini, we derive limits on the DM 
lifetime for several final states. Typically, those limits come 
from constraints on new sources of ionization and heating 
stemming from the products of DM interactions. Here we 
pay special attention to DM particle decays. Our findings 
rely on several standard assumptions namely: (i) the DM life¬ 
time is large than the age of the universe; (ii) the DM particle 
accounts for the DM cosmological abundance; (iii) the DM 
particle decays fully to SM particles; (iv) the energy fraction 
which the DM particle deposit into the intergalactic Medium 
is determined by the transfer functions provided in Ref. M- 
The rate at which a given DM particle decay induces heat¬ 
ing and ionization of the baryonic component of the IGM is 
proportional to 


B. Antiproton Data 

Since DM decays produce, in principle, matter and 
antimatter in equal amount, antiprotons are an interesting 
target for indirect DM searches, due to the relative rarity of 
antimatter produced in astrophysical processes. In this section 
we revisit the procedure to place bounds on the DM lifetime 
using antiprotons data from AMS-02 ll22ll . Antiproton data 
can be used to set stringent limits on the DM lifetime, since 
DM decays should at some level produce a sizable amount 
of antiprotons, even for leptonic final states through the 
inclusion of electroweak corrections (radiation of a gauge 
boson which decays hadronicallvl ll23l . 

With recent AMS-02 precise measurements of the antipro¬ 
ton/proton ip/p) fraction for energies up to 450 GeV one can 
derive new restrictive limits on the DM lifetime for several 
decay modes such as bb, WW and pp, since no evidence 
of new sources of antiprotons were found in the data[^ We 
follow Ref. Il25l . where constraints on the DM annihilation 
cross section were derived using an older data set for the 
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' See Ref. ES where a claim was put forth about excess antiprotons in the 
PAMELA data. 
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FIG. 1. Model independent bounds on the DM lifetime. Leftmost: CMB limits for the ee (gray), (green), rr (orange), bb (blue) and WW 
(magenta). Left/Center: Limits from AMS-02 using antiproton data for the bb,WW and jj,fi channels. Right/Center: Bounds on the spectral 
line emission from DM decays. Right: Constraints on the three-body f fv decays mode using measurements of neutrino detectors. 


Einasto DM profile based on the total antiproton flux. Here 
we will instead obtain limits on the DM lifetime with an 
NEW profile using the latest AMS data on the p/p ratio, and 
compare our finding with existing limits. 

The derivation of limits based on antiproton data is sub¬ 
ject to large astrophysical uncertainties associated with p pro¬ 
duction, propagation and solar modulation. Here, we employ 
the standard set of Min-Med-Max propagation models. Min- 
Med-Max represent values of diffusion parameters which pro¬ 
duce a mininum-to-maximum antiproton flux from a DM de¬ 
cay as shown in Table I (See Ref. Il2^ for a recent review). 
Current data seem to disfavor the Min propagation model ll27l 
and the Max-model induces arguably overestimated bounds, 
so we base our limits on the Med propagation model and an 
NEW DM distribution. To obtain limits on the DM lifetime 
we first solve the cosmic-ray transport equation in the Galaxy, 
in a steady state condition for the number density of antipro¬ 
tons (/p) per unit of kinetic energy T, 

Qp(T,rO + V-[iT(r,fOV/p-14(r)/p]-2M(z)r,„„/p = 0, 

(4) 

where we have neglected energy losses and re-acceleration 
processes. We describe below the physical meaning of each 
of those components: 

(i) The first term refers to the primary production of an¬ 
tiproton from DM decays expressed as. 


Qdm 


P p 
Mdm j ^ dK 


(5) 


where p is the DM halo profile assumed to be NEW, and 
dN/dK is the energy spectrum generated using PPPC4DM 

m. 


(ii) The second term accounts for the diffusion of cosmic- 
rays through their propagation in the interstellar medium. It 
is typically assumed to be constant in the diffusion zone and 
often parametrized in terms of the particle rigidity (momen¬ 
tum/atomic number) as follows. 


IC{K)=KoP{p/zY, (6) 

with P = v/c, where the normalization {Kq) and the spectral 
index {5) are associated with the properties of the interstellar 
medium and derived from measurements of the primary-to- 
secondary flux ratios of cosmic-rays such as Boron to Carbon 
||23, and obviously Z =1 for antiprotons. 

(iii) The third term refers to the convection mechanism 
which accounts for the drift of charged particles away from 
the disk, assumed to be infinitely thin with a half-height of 
lOOpc 1261 . induced by the Galactic Wind with a characteris¬ 
tic velocity Vconv and spatially constant in the diffusion zone, 
i.e., Vconv = sign{z)Vc. Departures from the thin disk as¬ 
sumption lead to one order of magnitude changes in the final 
limits as one can see in Pig. 6 of Ref. l28l . 

(iv) The fourth term represents the annihilations of antipro¬ 
tons with the interstellar gas which is proportional to 


Tpp — {jtH “F 4 ^ (7) 

where, (jpp = 0.661(1 -F 0.0115 

forT < 15.5 GeV and cTpp = 0.036T-°'56forr > 15.5 GeV 
ll29l 1^ . In Eq|7] we assumed that the helium-antiproton 
annihilation cross section is simply a rescaling of the proton- 
antiproton ||28]| . 

We now have all ingredients to solve the transport equation 
and to compute the astrophysical and DM decay predictions 
for the antiproton flux. A final physical effect, solar modula¬ 
tion, affects the prediction of the antiproton flux at the Earth’s 










atmosphere at energies below ~ 20 GeV, as result of the solar 
cosmic-ray wind and magnetic field. 
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We take into account this effect using the force-field ap¬ 
proximation, which determines the antiproton flux as a func¬ 
tion of the kinetic energy of the antiproton at the atmosphere 
(Tat) by re-scaling the interstellar flux which depends on the 
antiproton kinetic energy (T^s) as follows ED: 




( 8 ) 


with Tat = Ti, 
in Table I. 


6 ^, where (fTp is the Fisk potential as given 
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The combination of the bounds listed above are shown 
in Fig.l, third panel. Note that some operators discussed 
in our work give rise to Z 7 and /17 lines: for those we 
correct the energy of the gamma-ray line using the relation 
— nT'V/M jjm), where my is either the Z or 
Higgs mass, and divide the lifetime by a factor of two since 
we have one photon in the final state instead. 


D. Neutrino Data 


Using a data-driven model to account for the proton flux 
for the energy range of interest as presented by the PAMELA 
collaboration 021 . which is well fitted by a Fisk potential 
(f)^ = 0.7 (see fifth column of Table 1 of ESl l. one can Anally 
compute the total p/p ratio from primary and secondary 
production processes, as discussed above, and enforce the 
condition that the predicted p/p ratio does not exceed the 
ratio measured by AMS-02 data ESI at 95% C.L, for the 
specific choices of DM lifetime and mass; That results in 
constraints on the DM lifetime versus mass plane shown in 
the second panel of Fig. 1. Our results were obtained with 
PPPC4DMID code ESI . 

Notice that our limits are competitive with existing ones 
derived using PAMELA E^ and AMS-02 ESI . In particular, 
our limits are mildly similar to Ref ESI which performed a 
through analysis by including several energy loss processes 
we ignored. 


C. Line Searches 

If bright enough to be distinguishable from background 
continuum emission, and if morphologically diffuse E4l . 
gamma-ray spectral lines are known as a veritable smoking 
gun for DM annihilation or decay signals. Several experi¬ 
ments have searched for line emission at energies between 
10“^GeV up to 400 GeV: (i) Chandra and X-ray Multi-Mirror 
Mission (XMM) X-ray telescopes cover the 0.007 keV-12 
keV range ESl |36l; (ii) High Energy Astronomy Observa¬ 
tory (HEAO) accounts for the 3-48 keV E21; (iii) INTEr- 
national Gamma-Ray Astrophysics Laboratory (INTEGRAL) 
the 20k eV-7 Me V ESl ; (iv) The Imaging Compton Telescope 
(COMPTEL) along with The Energetic Gamma Ray Experi¬ 
ment Telescope (EGRET) screens the MeV-lOOGeV E9l : (v) 
Fermi-LAT covering energies up to 462 GeV ll40ll4Tll . Here, 
we simply combine all of those constraints. We point out that 
we make use here of the latest Fermi-LAT limits on the flux at 
95% C.L for the 180° region centered at the Galactic Center 
(R180), as described in Ref. Il40ll42l . Limits on the DM life¬ 
time are obtained after solving for the DM lifetime using the 
differential flux equation. 


Neutrino detectors are sensitive to DM decays and have 
been used to place limits on the DM lifetime. For a NFW 
profile for Pq — 0.4, the full sky differential neutrino flux 
from DM decays reads ll43l . 
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Another source of neutrinos from DM decays stems from 
cosmological decays of DM particle producing a diffuse neu¬ 
trino flux from decays at all redshifts which reads: 


d^u _ Ppm 1 j'°° dz 

dE^ AllM dm Tdm Jo ^ \ Hiz)) 

( 11 ) 

where pdm is the cosmological DM density, T[{z) = 

+ flm(l + z)^ is the expansion rate of the universe, 
and s{E^, z) is the universe opacity to neutrinos obtained in 
Ref il. The neutrino oscillation probabilities in vacuum 
is assumed to remain unchanged at the detector. Thus the 
primary neutrino flux from a specific flavor is redistributed 
equally into all neutrino flavors, so that the number of ex¬ 
pected events is given by. 


E, 


dN, 

' dE^ 




Nexp = (time x AH) ^ f J^j^Aeff{E^)dE^. 

( 12 ) 

By comparing with the 95% C.L limits on the number of 
events observed, constraints on the DM lifetime for two body 
decays were derived: (i) Ref. Il45ti47l used AMANDA and 
Super-K data; (ii) Ref. Il48l analyzed recent ICECUBE data; 
However, operator 014 in Table II induces three body decays 
if fv). Hence, we take the limits from those references and 
use PYTHIA 6.4 ll48l to derive the corresponding bounds on 
three body decay as shown in Pig|^ rightmost panel. We em¬ 
phasize that the 95% C.L. limits were obtained by requiring 
the theoretical predictions not to overshoot the data in any 
point. 
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Antriproton Propagation Model 
MED 5 = 0.7 Ko = O.OlUKpc^ /Myr 
Vc = 12km/s L = 4 Kpc = 0.7Gy 

TABLE I. Propagation model parameters: <5 and Kq are the spectral index and normalization that go into Eq.6; 14 is the wind velocity in Eq.7; 
L is the half-height of the cylinder with 20Kpc radius which is used to model the antiproton diffusion. We assumed the proton Fisk potential 
to be equal to the antiproton which is a good approximation as discussed in Ref. |28I . 
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FIG. 2. Model independent bounds on the DM lifetime. Left: Bounds on the spectral line emission from DM decays. Right: Constraints on 
the three-body ffv decays mode using measurements of neutrino detectors. 


E. Gamma-ray data 

Observations of the continuous emission of gamma rays 
give rise to stringent limits on the DM lifetime. Here we 
employ limits derived from: (i) the extragalactic gamma-ray 
background, as derived in Ref. Il49ll . which postulates that 
the sum of the isotropic component from blazars (making 
up nearly 70% of the total intensity), star-forming galaxies 
(SFGs), misaligned active galactic nuclei and DM decays not 
exceed the measured flux at 95% C.L (Fig.4 of Ref. Il49l k (ii) 
limits from Fermi-LAT observations of eight galaxy clusters 
at gamma-ray frequencies, in 10° x 10° squared regions cen¬ 
tered on the clusters EQl. We do not duplicate those results 
here, but we use them in what follows to derive our bounds. 
For several discussions related to gamma-ray constraints and 
other topics which provide complementary limits we point to 

RefsElEIMl- 


III. BOUNDS ON THE DARK MATTER LIFETIME 

As discussed in the previous section, we assume that 
continuous global symmetries are broken due to gravita¬ 
tional effects (see Ref m otherwise); in the presence of 
a continuous global symmetry, one should thus consider 
Planck-suppressed effective operators which break continu¬ 
ous global symmetries, leading to metastable DM particles. 
Limits on the lifetime of DM particles from observations in 
a broad range of frequencies thus allow us to derive general 


constraints on these operators in settings that advocate 
continuous global symmetries to stabilize DM candidates. 

We list in Table II a set of dimension-five gauge and 
non-gauge invariant operators that violate continuous global 
symmetries and induce DM decay. We point out that our 
list is not exhaustive or complete; rather, it should serve as 
a proof of principle, since it includes operators mediating 
several decay modes that produce significant continuum 
gamma-ray emission, spectral lines, antiproton, charged 
leptons and neutrino fluxes; in addition, the set we consider 
encompasses a variety of DM particle quantum numbers. 
Here we show results for the gauge-invariant operators only, 
while in the appendix we present results for the non-gauge 
invariant ones. Notice that other functional forms are pos¬ 
sible for Planck-scale operators violating continuous global 
symmetries, for example because of the effects of D-brane 
instantons, see e.g.Ref. CD. 

In the Table, we have introduced the dimensionless cou¬ 
plings Xi ~ 0{1), whose value depends on the unknown 
mechanism for the quantization of gravity. As we argue be¬ 
low, the precise values of Xi are irrelevant to our conclusions, 
but we keep the A^’s as free parameters and obtain our limits 
in the A vs DM mass plane. For each of the Planck-suppressed 
operator, we apply the most stringent limit on the DM lifetime 
for a given particle mass. Our results are collected in Figures 
3-4. In Figure 3 we show the result for gauge invariant oper¬ 
ators, whereas in figure 4 for non-gauge invariant one. Notice 
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Interaction Term 
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04 
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05 
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012 
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TABLE II. Dimension-five Planck-suppressed operators potentially 
inducing the decay of the DM particle in models with continuous 
global symmetries. 01-013 refer to scalar DM (S), 014 to fermion 
(ip), and 015-016 to vector (V) DM particles. In the table, r = ±1, 
H is the Standard Model Higgs, H the isospin transformation of H, 
Z and W are the neutral and charged weak gauge bosons, is the 
electromagnetic tensor, is the gluon tensor, and rep¬ 
resent the parts of the SU(2) tensor which involve the physical Z and 
W gauge bosons respectively, / and L represent SM fermions. The 
list includes both gauge invariant and non-gauge invariant operators 
that include a DM particle field singlet under the SM gauge group, 
but we will restrict our constraints to the gauge invariant ones. See 
however the appendix for results concerning the non-gauge invari¬ 
ant operators. The dark matter particles which come in multiplets 
higher dimensional operators are required giving rise to much larger 
lifetimes avoiding our limits. 


that several bounds are truncated at some DM mass due to the 
lack of data at lower energies. Moreover, the sudden change 
in the behavior of the limits in the figures has to do with the 
shift in the data set used to constrain a given effect operator. 
For instance, for operator Ol, for masses below 10 GeV we 
had to shift from gamma-ray to CMB data that yields much 
weaker constraints, accounting for the abrupt change in the 
limit at Mjjm = 10 GeV. 

Figures 3 shows that models that advocate the presence 
of continuous global symmetries to stabilize scalar DM 
candidates might produce a line emission with a very short 
lifetime (through operator 02/03), well below the age of the 
universe, thus ruling out DM masses larger than 100 keV. It is 
clear that for any scalar DM operator the whole electroweak 
WIMP range as well as the large mass range of warm DM is 
ruled out, since only for M^m ^ 100 keV are couplings of 
order one achieved. 


As for fermionic DM candidates stabilized by continuous 
global symmetries, operator 014 arises naturally at the Planck 
scale, yielding an appreciable neutrino and cosmic-ray flux as 
displayed in the right panel of Fig.3. The result clearly shows 
that we are able to exclude DM masses above 100 MeV. 
Lastly, in models where vector DM particles are stabilized via 
the existence of continuous global symmetries, operator 015 
would automatically be present at the Planck scale, leading to 
DM decay into fermion pairs. After employing a combination 
of the bounds shown in Fig.l, we And that masses larger than 
10 MeV induce cosmic-ray and gamma-ray fluxes that exceed 
the measured values. Conclusions regarding the remaining 
operators in Table II can be straightforwardly drawn. 

In summary, we And that dimension-five effective opera¬ 
tors at the Planck scale make continuous global symmetries 
problematic to stabilize DM particles outside very special, re¬ 
stricted mass ranges. In the next section we show how our 
bounds highly constrain several well-known models in the lit¬ 
erature. 


IV. CONCRETE MODELS 

In what follows, we work under four assumptions: (i) 
global symmetries are in fact broken at the Planck scale;(ii) 
0 (1) non-renormalizable couplings are a good approximation; 
(iii) dark matter held is a singlet held; (iv) the dark matter den¬ 
sity distribution is well described by a NFW profile. 

It is important to keep in mind as caveats that despite the 
existence of plausible arguments ruling out global symme¬ 
tries at the Planck Scale, the validity of such arguments is still 
under debate. Moreover, in general continuous global sym¬ 
metries can be replaced by other symmetries such as discrete 
symmetries, circumventing our constraints. We assume 0(1) 
non-renormalizable couplings for the non-renormalizable op¬ 
erators mediating decay. Any departure from this assumption 
is easily factoring in by rescaling using the expressions for the 
width rate given in Table 3. 

We apply here the results from the previous sections to con¬ 
crete models, widely discussed in the literature where the dark 
matter held is a singlet held, for which our bounds are directly 
applicable. 

Left-Right Model 

We consider the left-right mirror symmetric model with the 
continuous global symmetry (7(1) b-l ® (7(l)xofRef 173] 
ED. There, the continuous global symmetry prohibits the term 
LHipn, where '0 is a fermionic DM candidate. However, 
such symmetry is generically violated at the Planck scale and 
therefore the 014 operator ought to exist, thus ruling out DM 
masses above 100 MeV. Ref. 1231211 also invokes the case 
of WIMP scalar DM protected by the continuous global sym¬ 
metry, but once again, as we see in the left panel of Fig.2, 
the entire corresponding WIMP mass range is excluded. As 
a result, the model described in Ref. CSEll does not appear 
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FIG. 3. Limits on the A as a function of the DM plane by enforcing that energy injected at dark ages does not distort the CMB power spectrum 
at 95% C.L., and the flux of neutrinos, cosmic-rays and gamma-rays induced by DM decays of each individual operator does not exceed the 
measured values at 95% C.L. The limits are based on the results presented in Figures 1-2. These bounds are applicable to any scalar, fermion 
or vector DM particle stabilized by a continuous global symmetry. From our results one may exclude: (i) scalar DM with Mdm > 100 keV; 
(ii) fermion DM with Mdm > 100 MeV; (iii) Vector DM with Mdm > 10 MeV. 


to have a plausible DM candidate. Unless the invoked con¬ 
tinuous global symmetry can be replaced by a Zn discrete 
symmetry of some sort, the model is strongly disfavored by 
data. 

Two Higgs Doublet Model 

In the original two Higgs doublet model no DM candidate 
is present. Nevertheless, if the second Higgs doublet is odd 
under a symmetry the CP-even scalar of that doublet can 
be a DM candidate. This is the case in the so-called Inert Two 
Higgs Doublet Model (I2HDM) ll75l . Recently, a continuous 
global symmetry has been proposed to replace the sym¬ 
metry ll76l . The authors focus on the spontaneous symmetry 
breaking of the continuous global symmetry through the vev 
of the second Higgs doublet and comment on the possibility of 
having an unbroken continuous global symmetry. Despite the 
interesting Higgs physics implications produced by the use of 
the continuous global symmetry m, operator 02 should be 
present at the Planck scale; thus, from Fig. 2 we conclude that 
DM masses larger than 100 keV are problematic along with 
the possibility of having viable WIMP DM candidates in the 
model. 

Singlet Fermion Model 

The minimal fermionic DM model studied in Ref. llTTlI ad¬ 
vocates a continuous global symmetry responsible for stabi¬ 
lizing a singlet fermion which yields the desired thermal relic 
abundance and is consistent with direct searches. The Planck 
suppressed effective operator 014, however, rales out the en¬ 
tire WIMP mass range. As we mentioned before, in princi¬ 
ple one could replace the continuous global symmetry by a 
discrete symmetry, since Zm is a subgroup of U{1). How¬ 
ever, the necessary discrete symmetry might turn out to imply 
a rather large and unnatural tuning of the model. 


Radiative See-Saw Model 

A radiative lepton model in which the charged lepton 
masses are generated at one-loop level whereas and the neu¬ 
trino masses at two-loop level has been proposed in Ref. ED. 
In this model the continuous global and Z 2 symmetry have 
been invoked and two DM candidates postulated. A singlet 
fermion, referred to as n' in Table I of Ref. Ga, is not odd 
under the Z 2 symmetry, and claimed to be a WIMP due to the 
presence of a continuous global symmetry. Similarly to the 
previous model, Planck-suppressed dimension-five operators 
discard such possibility. 

Zee-Babu Model 

The Zee-Babu model adds to the SM a singly-charged and a 
doubly-charged scalar ll7^ . Recently, an extension of the Zee- 
Babu model has been put forth by adding a singlet fermion 
which is stabilized by a continuous global U{1)b-l symme¬ 
try. This continuous global symmetry also forbids terms like 
LHN. There, the neutral fermion does not carry a lepton 
number so it is purely a neutral fermion. Nevertheless, as we 
discussed, this continuous global symmetry does not hold up 
to the Planck scale and consequently the operator 014 arises, 
inducing an excess production of neutrinos, gamma-ray and 
comic-rays, which results into the exclusion of DM masses 
below 100 MeV, in tension with what presented in Figs. 1-2 of 
Ref. f79l . 

3-3-1 Models 

3-3-1 models refer to «S'(3)c <8 ) 5'(3)l 8 U(l) 7 v gauge ex¬ 
tensions of the Standard Model lISOl . In Ref. lISTl a contin¬ 
uous global symmetry with the purpose of avoiding undesir¬ 
able mixing among the gauge bosons and of guaranteeing that 
the lightest particle charged under the continuous global sym¬ 
metry be stable. Both a complex scalar and a heavy Dirac 
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fermion were studied as potential WIMP DM candidates. In a 
similar vein to what discussed above, the WIMP mass regime 
in this model is in jeopardy due to the aforementioned gravity 
effects. In the urge of preventing the use of continuous global 
symmetries in the model. Ref. Il82l proposed adding an extra 
gauge symmetry, which would completely change the associ¬ 
ated DM phenomenology. 

V. CONCLUSIONS 

Based on general lessons from quantum mechanics and 
general relativity, it reasonable to assume that no global sym¬ 
metries are allowed in a consistent theory of quantum gravity. 
We have presnted, in a model-independent approach, robust 
gamma-. X-ray, CMB, and cosmic-ray constraints on decay¬ 
ing DM particles, using a large set of data, including data from 
Fermi-LAT, AMS-02, Super-Kamiokande, Planck, WMAP9, 
AMANDA, and Icecube among others. We have then ap¬ 
plied those bounds to scalar, vector and fermion DM particles 
decaying through dimension-five Planck-suppressed effective 
operators. We stress that our findings are based on the follow¬ 
ing assumptions: 

(i) global symmetries are broken at the Planck scale; 

(ii) the operators mediating dark matter decay have 0(1) 
non-renormalizable couplings; 

(iii) the dark matter particle is represented by a singlet field; 

(iv) the dark matter density distribution in the Galaxy is 
well described by a NFW profile. 

Under these assumptions, we have derived the following 
constraints on the possible mass range: 

(i) scalar DM: Mdm ^ 100 keV ; 

(ii) fermionic DM: Mdm ^ 100 MeV ; 

(iii) Vector DM: Mdm < 10 MeV . 

Lastly, we have applied our limits to instances of parameter 
space of Left-Right, Two-Higgs Doublet, Singlet Fermionic, 
Zee-Babu, 3-3-1 and Radiative See-Saw models to conclude 
that the occurrence of DM particles in such models is generi- 
cally problematic outside the DM particle mass ranges listed 
above. Our results basically rule out the entire WIMP mass 
range for models invoking continuous global symmetries to 
stabilize the DM particles. 

We emphasize that our results rely on the critical assump¬ 
tions listed above, especially on the hypothesis that continu¬ 
ous global symmetries are violated at the Planck scale. Notice 
that if one uses non-renormalizable couplings smaller than 
unit, our limits can be simply rescaled using Table 3 con¬ 
taining the decay rates. Moreover, if one adopted a different 
density profile that would quantitatively change the accept¬ 
able DM mass ranges, but leave the overall conclusions un¬ 
changed. However, if the DM particle belongs to a multiplet, 
and thus only higher dimension operators induce the dark mat¬ 


ter decay, our limits would not be relevant. In summary, our 
bounds and conclusions can be subject to significant changes. 

Mambini, Profumo, Queiroz 2015 



FIG. 4. Limits from non-gauge invariant operators on A as a function 
of the DM mass, obtained by enforcing the 95% C.L. bounds from 
the various observational probes discussed above. 
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VI. APPENDIXrDECAY WIDTHS AND NON-GAUGE 
INVARIANT OPERATORS 

In Table III we list the decay width associated with each 
operator discussed in the manuscript; in addition we discuss 
some non-gauge invariant operators (06-09-012-013) that 
might appear in more complex setups such as non-Abelian 
theories lf83l . We also show the limits stemming from such 
operators in Fig|^ 
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